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Abstract 



We study a version of compact directed percolation (CDP) in one 
dimension in which occupation of a site for the first time requires that 
a "mine" or antiparticle be eliminated. This process is analogous to 
the variant of directed percolation with a long-time memory, proposed 
by Grassberger, Chate and Rousseau [Phys. Rev. E 55, 2488 (1997)] 
in order to understand spreading at a critical point involving an infinite 
number of absorbing configurations. The problem is equivalent to that 
of a pair of random walkers in the presence of movable partial reflectors. 
The walkers, which are unbiased, start one lattice spacing apart, and 
annihilate on their first contact. Each time one of the walkers tries to 
visit a new site, it is reflected (with probability r) back to its previous 
position, while the reflector is simultaneously pushed one step away 
from the walker. Iteration of the discrete-time evolution equation for 
the probability distribution yields the survival probability S{t). We 
find that S{t) ~ t"*^, with 5 varying continuously between 1/2 and 
1.160 as the reflection probability varies between and 1. 
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I. INTRODUCTION 



Models that can become trapped in one of an infinite number of absorbing con- 
figurations (INAC) exhibit unusual spreading dynamics at their critical point. The 
most intensively studied model of this kind is the pair contact process (PCP) 
INAC appears to be particularly relevant to the transition to spatio-temporal chaos, 
as shown in a recent study of a coupled-map lattice with 'laminar' and 'turbulent' 
states, which revealed continuously variable spreading exponents [0]. 

Anomalies in critical spreading for INAC (such as continuously variable criti- 
cal exponents) have been traced to a long memory in the dynamics of the order 
parameter, p, arising from a coupling to an auxiliary field that remains frozen in 
regions where p = @J§]. Grassberger, Chate and Rousseau (GCR) proposed 
that spreading in models with INAC could be understood more easily by studying a 
model with a unique absorbing configuration, but with a long memory of its initial 
preparation. 

The GCR model is a variant of bond directed percolation (DP) in which bonds 
connecting to "virgin" sites (i.e., that have never been occupied), have a transmission 
probability, q, that may differ from the value, p, for bonds to "used" sites. Used 
sites follow the usual DP rule. If site x has been occupied previously, then the 
probability that x is occupied at time t + 1 is p if either x — 1 or x + 1 (but not both) 
are occupied at time t, p{2—p) if both sites are occupied at time t, and zero if neither 
is occupied. For virgin sites the parameter p is replaced by q. The dynamics begins 
(as in all spreading experiments) with activity restricted to a small region of the 
lattice. Grassberger et al. found in simulations that in 1-1-1 dimension, the critical 
point remains at Pc = 0.644701, the standard bond DP value 0, independent of q. 
They concluded that for q < pc, the survival probability S(t) decays faster than any 
power of t, at the critical point p=Pc- 

In this work we study one-dimensional compact directed percolation (CDP) [§0, 
so called because gaps cannot arise within a string of occupied sites. Being exactly 
soluble, CDP provides a valuable test for ideas on scaling in absorbing-state phase 
transitions. For example. Odor and Menyhard recently found a continuously- variable 
survival exponent for CDP confined to a fixed parabolic region |]TD[. 

The rules of standard CDP are as for DP, described above, except that if both 
X — 1 and X + 1 are occupied at time t, then x must be occupied at time t + 1. 
(Note that CDP possesses two absorbing states: all vacant, and all occupied.) If the 
process starts with only a single occupied site, the state at any later time is specified 
by the positions of a pair of random walkers, wi and W2, which mark the extent of 
the occupied region. (Specifically, the occupied sites are: wi + l,wi + 2, ...,W2-) If 
we take the origin as the position of the original "seed" particle, then wi = and 
W2 = 1 at t = 0. The stochastic evolution of W2 is given by 

W2{t+1} = < },{ ^ ^ (1) 

yw2{t) — i w.p. 1—p ^ ' 

while for w\ the roles of p and \ —p are interchanged. Thus the length Yii) = 
W2{t) — Wi{t) of the occupied region itself executes a random walk with transition 
probabilities: 
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w.p. (1— p)^ 

Y{t + l) = {Y{t) w.p. 2p{l-p) (2) 




w.p. p"^ 

The state Y = (all sites vacant) is absorbing. In this work we focus on the case 
p = l/2 since for smaller (larger) values Y{t) is attracted to (driven away from) the 
origin. 

Well known results on random walks fill imply that for (j9= 1/2), the survival 
probability S{t) ~ t^^^^, while the mean-square displacement, if the walker has not 
hit the origin up to time t, follows {Y'^{t))s ~ t. (The average is over trials that 
survive until time t or longer.) The latter implies an active region (in surviving 
trials) of extent ~ t^^'^, so that the mean number of occupied sites, averaged over all 
trials, is n ~ t^. In the usual notation of absorbing-state phase transitions |]T2|,|13[, 
these results imply the exponent values 6 = 1/2, z = 1 and rj = 0. (The exponents 
are defined via the relations P ~ t"^, R"^ ^ t^, and n ^ .) These values satisfy the 



expected hyperscaling relation for a compact growth process in d dimensions 



s + v = Y- (3) 

Now we introduce a memory effect in CDP along the lines proposed in Ref. [Q. 
Suppose that initially sites other than the origin harbor static "antiparticles" (or 
"mines"), independently with probability r. If site x has a mine, then the first 
particle to venture there is destroyed, and along with it the mine, so that in future, 
site X can be occupied as in normal CDP. In terms of the random walkers Wi and 
W2, a mine is effectively a reflecting boundary: the first time Wi attempts to visit 
site X (mined), it is reflected back to x + 1, and at the same time the reflector moves 
to X — 1; similarly, W2 will be reflected back to x — 1 on its first visit to x, if it 
harbors a mine. For r > our model represents the spread of activity into a hostile 
environment, for example the advance of a bacterial colony in a medium, with a 
preliminary contact facilitating expansion into new regions, or, similarly, the spread 
of a political viewpoint in an initially skeptical population. 

We recently studied a simplified version of this problem, involving a single ran- 
dom walker on the nonnegative integers |T3. The walker is unbiased, and starts at 



X = 1, with a; = absorbing. The reflector is initially at x = 2. On each visit to 
a new site, the walker is reflected with probability r, and the reflector moves for- 
ward by one site. Asymptotic analysis of the probability generating function shows 
that the survival probability exponent varies continuously with r: S = (l+r)/2. 
In this work we analyze the two-walker problem defined above, corresponding to a 
spreading CDP process. 

The remainder of this paper is organized as follows. In the following Section we 
show how CDP with reflectors can be represented (despite the long memory) as a 
discrete-time Markov process. We proceed to define an appropriate state space and 
the associated transition probabilities. In Sec. Ill we analyze the results of numerical 
iteration of the probability evolution equations, yielding precision estimates of the 
critical exponent 6 and other asymptotic scaling properties. Sec. IV presents a 
summary and discussion. 
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II. MODEL 



To investigate the scaling properties of CDP with reflectors, it is convenient to 
enlarge the state space to include the positions of the reflectors; this renders the 
process Markovian. [The process {wi{t),W2{t)) is evidently non-Markovian.] We 
consider CDP in one dimension, starting from a single active site. The evolution of 
the active region is represented by the motion of a pair of unbiased random walkers, 
wi and W2. Each time a walker tries to jump to a new site it pushes a reflector 
or R2) to the left or right, and in the process the walker is reflected back with 
probability r; it remains at the new position with probability f = 1— r. The generic 
conflguration is: 

Ri wi . . . W2 R2 

Due to translation invariance, we require only three variables, x, y and z, deflned 
as follows: 

X — wi — Ri — 1, 
y ^ R2 - W2 - 1, 
z = R2- Ri-2. 

Then the distance between the walkers (i.e., the number of occupied sites), is z—x—y; 
x + y = z is the absorbing state. We start with the walkers a unit distance apart, 
and the reflectors one lattice spacing away from the walkers, so that, initially, z = 1 
and X — y — 0. At each time step the walkers jump to the left or right with equal 
probabilities, z is nondecr easing, with 1 < z < 2t + 1, since the separation between 
reflectors can increase by at most two spacings at each step. 

Let P{x,y, z;t) denote the probability of state {x,y,z) at time t. Transitions 
{x, y, z) — > {x', y', z') may be grouped into three classes. The simplest is for x and y 
both greater than zero. Then z cannot change, since the walkers do not encounter 
the reflectors, and we have {x, y, z) — > (x', y', z) with x' — x^l and y' — y±l; each 
of these has a transition probability W = 1/4. Next consider y > x = 0. There 
are six possible transitions, listed, along with their probabilitites, in Table 1. (The 
transition probabilities for x > y are obtained by noting that W is symmetric under 
the simultaneous interchange of x and y and x' and y'.) 
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x' 


y' 


z' 


W 







z+l 


r/4 


1 




z+l 


r/4 





y-i 


z+l 


r/4 


1 




z+l 


r/4 


1 


y+l 


z 


1/4 


1 


y-1 


z 


1/4 



Table I. Transition probabilitites for ^/ > x = 0. 
Finally, for x = y = 0, there are eight possible transitions, as listed in Table II. 
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z+l 
z+l 
z+l 
z+2 
z+2 
z+2 
z+2 
z 



W 



r/A 

f/4 

r/2 

fV4 

rr/4 

rr/4 

rV4 

1/4 



Table 11. Transition probabilitites lox y — x — ^. 

(Note that there are two distinct routes to the state (1, l^z + l): both walkers may 
jump to the left, with wi reflected back, or both may jump to the right, with W2 being 
reflected; each of these events has a probability of r/4.) Any move yielding x' + y' > 
z' represents a transition into the absorbing state. Starting from P(x, y,2;;0) = 
(^^,i^x,o^?/,0) we can iterate the above transition probabilities to find P{x,y, z;t). 

In the three- variable representation, the evolution is confined to an infinite wedge 
bounded by the planes x = 0, y = 0, and x+y = z. The latter plane is absorbing, while 
the first two allow upward transitions (from z to z + l or z + 2). Between vertical 
transitions, the process is confined to the triangle x > 0, y > 0, x+y < z; away from 
the boundaries, the evolution is that of an unbiased lattice walk with steps between 
second-neighbor sites, on Z^. 

Suppose the process has just entered a given plane of constant z from below. 
Its continued survival is equivalent to the event that it touches either the x or the 
y axis, and makes a further vertical transition, before touching the line x+y — z. 
Thus survival of the process is related to the splitting probabilities for exiting a two- 
dimensional triangular region via the different edges. (Note that the x and y axes 
are partly refiecting.) 

The above transition probabihties define what we shall refer to as the "two-step" 
model, in which both walkers jump at each time step, in correspondence with the 
original CDP problem. One may define a simpler "one-step" model, in which only 
one of the walkers (chosen at random, with equal likelihood) jumps at each step. 
Since the set of transitions is somewhat reduced (there are no transitions from z to 
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for example), this version woTild appear to be more amenable to analysis. We 
expect the two versions to have identical asymptotic scaling properties. 

The transition probabilities for the one-step process are as follows. If x and y are 
both greater than zero, there are four possible transitions, (x, z) ^ {x ± 1, y, z), 
and (x, y,z) — > {x,y ±1, z), each with probability W — l/A. For y > x — 0, the five 
possible transitions are listed in Table III. 
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Table III. One-step model: transition probabilitites for y > x = 0. 

Finally the transitions for the case x = y = Q for the one-step model are given in 
Table IV. 
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Table IV. One-step model: transition probabilitites lor x — y — Q. 

The evolution of the one-step process is again confined to the wedge described 
above, and (in each plane) to the same triangular region as the two-step process. The 
principal differences are that, away from the boundaries, the process corresponds to 
a simple random walk with jumps between nearest neighbors, and that all vertical 
transitions are from z to z+l. 

The probability distribution evolves via 

P{x,y,z-t+l)^ J2 W{x,y,z\x',y',z')P{x',y',z';t). 

x',y',z' 

Results from iteration of this equation are discussed in Sec. III. 

A. Calculational scheme 

Given the symmetry of the transition probabilities under exchange of x and y 
(and, simultaneously, of x' and y'), it follows that if we start from a symmetric 
distribution, P{x,y,z) = P{y,x,z), as is the case here, then this property will be 
maintained throughout the evolution. This allows us to reduce the number of states 
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by roughly half: wc need only study x < y. The presence of the absorbing state 
implies that y < z~l, and, therefore, < x < min[y, z — y — 1]. 

Since states with x > y are not considered explicitly, we must modify the iteration 
of the evolution equation as follows: 

(1) The contribution to P{y', y', z'; due to a transition (,x, y, z) (y', y', z) with 
x<y should be doubled, to take into account the corresponding contribution due to 
{y,x,z) — > {y',y',z), which is not represented explicitly in the dynamics. Similarly, 
a transition from {x, y, z) to the absorbing state should have its weight doubled, if 
X <y. 

(2) In a transition {x,y,z) — > {x',y',z'), with x < y and x' > y', the contribution 
should instead be added to P{y',x',z';t + 1), to include the mirror process, which, 
again, is not represented explicitly. 

These rules are summarized in Table V, which gives the weights associated with 
each transition, given that states with x > y are not represented explicitly. 



From x = y to 


weight 


x' = y' 


1 


x' < y' 


1 


x' > y' 





absorbing 


1 


From x<y to 


weight 


x' = y' 


2 


x' < y' 


1 


x' >y' 


1 for {y', x') 


absorbing 


2 



Table V. Transition weights. 



The zero entry ior x— y^ x' >y' means that such transitions are ignored. 

III. RESULTS 

We have iterated the discrete-time evolution equation derived above numerically. 
To iterate the two-step process for tm = 2000 — 5000 time steps, one requires values 
of z of up to 200 - 320, and of x and y up to 110 - 230, depending on r. (The larger 
r is, the less rapidly the process spreads, and the smaller the arrays need be. The 
required size scales, naturally, as y/t^. For t„j = 2000 the iteration requires about 
20 min. to 1 hour of cpu time on an alpha workstation.) For the one-step process we 
use an upper limit of 250 for all three variables, which proves more than sufficient for 
tm — 2000. The CDP plus reflectors problem is, of course, easily studied via Monte 
Carlo simulation. But we have found that numerical iteration furnishes an order of 
magnitude higher precision than direct simulation, for the same expenditure of cpu 
time. 
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For each r we calculate the survival probability S{t), the first and second mo- 
ments, {Y)t and {Y'^)t of the extent of the active region (i.e., the distance between 
the walkers, z—x—y), and the probabihty distribution P(Y) at tm- (The distribution 
and moments of Y are taken over the surviving sample at time t. With the array 
sizes mentioned above we determine the quantities of interest to a precision of better 
than one part in 10®.) 

S{t) and the moments of Y are found to follow power laws, 

S{t) ~ (4) 



{Y)t - t^' (5) 

{Y')t ~ t"'^ (6) 

(Here the subscript 's' denotes surviving sample; the exponent rj = rjs — 6.) 

Precise estimates of the exponents are obtained by studying local slopes, for 
example 6{t) = din P/ dint, and similarly for the other exponents. Based on experi- 
ence with the random walk with movable partial reflectors |T3, we expect a generic 
correction to scaling exponent of 1/2; we therefore plot the local slopes versus 
Such plots (see Fig. 1), are roughly linear, but show a certain degree of curvature, 
indicating (as is to be expected) that corrections of order are still significant. 
The local slope data are fit nearly perfectly by a quadratic form in t~^/^; the inter- 
cept yields our estimate for the critical exponent. An exception is the case r = 1, 
for which the correction to scaling exponent appears to be 1. (This may be seen 
explicitly in the case of a single random walk with a partial movable reflector ||15|| .) 
For r = 1 we derive our estimate for 6 from an analysis of the local slope as a function 
of t~^. The extrapolated values for 6 are very stable under changes in the interval 
used (e.g., r^/^ < g.l, or r^/^ <- o.04), and in (2000 or 5000 time steps). We 
estimate the uncertainty of extrapolation as < 2 x 10~^. This is supported by our 
results for the one-step process: the estimates for 6 differ from those for the two-step 
model by at most 0.0005. (For r = 1 for example, we find 6 = 0.1597, compared 
with 6 = 0.1595 in the two-step case.) 

The analysis described above yields 77^ = 1/2 to within one part in 5000. Thus 
the only independent exponent is S. Our results for 6{r) are given in Table IV. As 
shown in Fig. 2, S appears to vary linearly with r. A simple linear expression, 

S='-^'A (7) 

2 3 ^ ^ 

reproduces the data to within 6 parts in 1000. The simplicity of this expression, 
and its similarity to the single- walker result, 6 = l/2 + r/2, lead us to adopt Eq. 
(0) as a conjectured exact formula. Analysis using least-squares fitting suggests, 
however, that 6{r) is weakly nonlinear. We obtain an excellent fit to our data using 
an expression of the form, 

6 =^ + ar + br^, (8) 
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with a = 0.6434 and b = 0.0167. The typical errors associated with a purely 
linear fit are on the order of 2 x 10~^, about an order of magnitude larger than the 
uncertainty in S. (The typical error for the best quadratic fit is about 4 x 10~^.) 
The nonlinear dependence of 5 on r therefore appears to be real, not just an effect 
of finite numerical precision. It is conceivable, nevertheless, that corrections to 
scaling introduce systematic errors in the numerical analysis, giving rise to apparent 
nonlinearities. We defer the verification of Eq. (|^, which would appear to require 
either an analytic solution or improved numerics, to future work. 
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0.9591 
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2174 
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0.8 


1.0258 
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1722 
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2411 


0.9 


1.0927 
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1316 
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2386 
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1.1595 


1 


0947 


1 


2362 



Table VI. Numerical results. 



Table VI also contains results for the amplitude A of the mean activity, defined 

via 

(F) ~ At^/\ (9) 
and for the asymptotic moment ratio 

The amplitude decreases smoothly with r as shown in Fig. 3. m is a measure 
of the shape of the position distribution. For r = our numerical estimate is 
consistent with A/tc = 1.27324..., as expected for Brownian motion on the line 
with the origin reflecting, for which the asymptotic probability density is Py(x) = 
(x/a^) exp[— x^/2cr^], with = t. The ratio decreases steadily with r, but not 
by very much (see Fig. 3), showing that the random walk result still serves as a 
reasonable approximation for r > 0. The moment ratio m takes the same values in 
both versions of the process, confirming that it is a universal quantity. 

The effect of the refiectors is clearly evident in the distribution in the number of 
active sites, P{Y). Fig. 4 compares P{Y) for r = 0, 0.5, and 1, (for t = 2000 time 
steps, two-step process), showing that the distribution shifts to smaller Y values 
with increasing refiection probability r. (Note that for r = 0, P{Y,t) = for Y + t 
odd. We have therefore multiplied the distribution for r = by one half, to facilitate 
comparison with the other cases.) Despite the changes in form, the tail of the 
distribution remains Gaussian in all cases. 
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IV. DISCUSSION 



We have studied spreading in compact directed percolation on a one-dimensional 
lattice at its critical point, modified so that activation of virgin sites is less probable 
than reactivation of a previously active site. The problem is equivalent to a pair of 
random walkers subject to movable partial reflectors. Two variants are considered: 
the two-step process, in which both walkers move at each step, and a one-step 
process, in which, at each step, only one walker (selected at random) jumps. We 
study these processes via exact numerical iteration of the probability distribution 
for finite times (< 5000 time steps). 

We find that the survival probability critical exponent 5 varies continuously with 
r. Our results indicate a weak nonlinearity in the function S{r), despite the fact that 
the data are rather well represented by a simple linear expression, Eq. Since, 
in the case of a single random walker subject to a partial movable reflector [|T3|, we 
found a strictly linear dependence of the survival probability exponent on r, this 
nonlinearity is somewhat surprising. On the other hand, the present problem is 
related to splitting probabilities on a two-dimensional domain (rather than on the 
line, as is the case for a single walker), allowing for a more complicated functional 
dependence. 

Our finding of a continuously-variable survival probability exponent is fully con- 
sistent with previous results for the single walker. One may, moreover, understand 
the fact that 6 increases more rapidly with r than for the single walker, since in GDP 
the spreading process feels the effects of two reflectors. In the absence of an exact 
analysis or rigorous argument, however, we have no qualitative understanding of the 
values for 6{r) that we have found numerically. This question, and the verification 
of Eq. d^), remain as interesting challenges for future work. 
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Figure Captions 

FIG. 1. Inset: survival probability S{t) for the two-step process, r = 0.8. Main 
graph: local slope 5{t) versus t'^^"^ for the same system, 'x' on y-axis: (upper) 
extrapolated value {6 = 1.0256); (lower) 6 = 1.0333 predicted by Eq. (|^. 

FIG. 2. Survival exponent 6{r) from iteration of probability distribution (points). 
The solid line is a quadratic least-squares best-fit to the data. Inset: best-estimate 
for 6 less the value predicted by Eq. (|^). 

FIG. 3. Upper panel: amplitude A of the mean activity as a function of r; lower: 
moment ratio m of the activity distribution. 

FIG. 4. Probability distribution P{Y) of the activity in the two-step model (con- 
ditioned on survival) after 2000 time steps, for (left to right) r = 1, 0.5, and 0. 
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